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Abstract Magnetic fields are fundamental to the dynamics of both accretion disks 
and the jets that they often drive. We review the basic physics of these phenomena, 
the past and current efforts to model them numerically with an emphasis on the jet- 
disk connection, and the observational constraints on the role of magnetic fields in 
the jets of active galaxies on all scales. 
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1 Introduction 

The topic of this workshop was the diverse astrophysical phenomena enabled by the 
strongest magnetic fields in the Universe. Gamma ray bursts and magnetars provide 
two examples of extreme phenomena where extraordinarily strong fields are manifest. 
However, the concept of “strong” is relative. A magnetic field may be weak in the 
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sense that the energy density associated with the field is small compared to other 
energies, e.g., thermal, rotational, gravitational. Such “weak” magnetic fields can, 
however, enable phenomena that otherwise would not be possible. In this contribution 
we consider just such an example, namely accretion disks and the jets they produce. 
Magnetic fields are fundamental to both accretion and outflows, and while in some 
specific cases the magnetic fields may be strong in the conventional sense of total 
energetics, for the most part the fields are likely to be relatively weak by those same 
conventional metrics. Nevertheless, magnetic fields play the central dynamic role. 

Gravity is the dominant force in creating the structure in the universe. When grav¬ 
ity is combined with angular momentum, the result is a disk. Since some angular mo¬ 
mentum will inevitably be present, disks are arguably the most basic type of object 
in the universe. Other structures, such as stars, result from disks that have somehow 
shed most of their angular momentum. Accretion is also a major astrophysical power 
source. As Lynden-Bell observed fLynden-Bell| 1969 ), the gravitational energy avail¬ 
able in black hole accretion greatly exceeds that available through nuclear reactions, 
with efficiencies anywhere from 6^-0% me 2 , versus less than 1% for nuclear fusion. 

Direct observational evidence for magnetic fields in accretion disks is very limited 
(measurements exist, e.g., for FU Orionis, Donati et al. |2005] ). The strongest case for 
magnetic fields is theoretical. We know from observations that disks accrete at rates 
that require considerable internal stress to transport angular momentum. As we shall 
discuss, this requires magnetic fields. 

Indirect evidence for magnetic fields in accreting systems comes from observa¬ 
tions of highly collimated outflows, or jets. Jets are one of the most striking signa¬ 
tures resulting from accretion, and are seemingly ubiquitous; they are associated with 
young stars, micro-quasars or X-ray binaries (MQs, XRBs), and active galactic nu¬ 
clei (AGN). The current understanding of jet formation is that outflows are launched 
by magnetohydrodynamic (MHD) processes in the close vicinity of the central object 


(Blandford and Payne| 1 982]|Pudritz and Norman|1983)|Pudritz et al.|2007)|Shang~et| 


al.|2007 1 It is furthermore clear that accretion and ejection are related to each other. 

Jets may be a consequence of accretion, but accretion may, in turn, be a consequence 
of a jet. Accretion occurs when gas loses some of its angular momentum, and one 
very efficient way to remove angular momentum from a disk is to eject it vertically 
into a jet or outflow ( |Wardle and Konigl|1993| |Ferreira and Pelletier|1993[ |Li| 1995 1 
Ferreira|1997| l. Not all details of the physical processes involved are completely un¬ 
derstood yet. 

In stellar sources, a central stellar magnetic field is surrounded by a disk carrying 
its own magnetic flux. Such a geometrical setup can be found in young stars, cata¬ 
clysmic variables, high-mass and low-mass X-ray binaries, and other micro-quasar 
systems. Systems such as AGN or ji -quasars have a black hole at their center, and a 
black hole cannot have an intrinsic magnetic field. The surrounding accretion disk, 
however, can support the magnetic flux needed for jet launching. In this case, the 
interaction of the black hole itself with the ambient field may launch a highly mag¬ 
netized Poynting flux-driven axial flow (Blandford and Znajek|1977 1. 

Although jets are found in a wide range of systems, from protostars to AGN, 
from the observational point of view protostellar and relativistic jets are quite differ¬ 
ent. For protostars we can measure Doppler-shifted forbidden line emission to obtain 
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gas velocities, densities or temperatures, and these provide a good estimate on the 
mass flux of the outflow. Evidence concerning the magnetic field, however, is mainly 
indirect: magnetic flares of large-scale protostellar stellar magnetic fields indicate 
that protostellar jets originating from a area close to the star, originate in an magne¬ 
tized environment. As noted above, in a few objects there are direct measurements 
of the strength or orientation of magnetic fields in the disks (e.g. Donati et al. 2005 


Steph ens et al.|2014| ). Concerning the jet magnetic field, direct evidence is less clear. 
So far an estimate is available for only one source, HH80/81, obtained from observed 
synchrotron emission ( Carrasco-Gonzalez et al.||2010[ l. The situation is completely 
reversed for relativistic jets. While there is plenty of information concerning the jet 
magnetic field structure, and some constraints on field strengths, speeds are harder 
to measure, there is no direct measure of the mass fluxes involved, and the magnetic 
properties of the accretion flow are not directly observable. We do not even know for 
certain whether the matter content of these jets is leptonic or hadronic (Section[5]i. 

While protostars and AGN are the most common jet sources known, the typical 
characteristics of jet systems — the existence of an accretion disk and a strong mag¬ 
netic field — are also present in other astrophysical sources, namely in cataclysmic 
variables, high mass or low mass X-ray binaries, and pulsars. Resolved Chandra 
imaging of the Crab ( Hester et al.|2002) and the Vela pulsar ( |Pavlov et al.|2003] l do 
indeed show elongated, highly time-variable but persistent structures emerging along 
the rotational axis of these systems. Whether these structures are similar to the jets 
observed in protostars and AGN or are features intrinsic to the pulsar wind nebula is 
not yet clear. Persistent jets are not widely observed in cataclysmic variables (either 
in the highly transient dwarf novae or the more stable nova-like objects), although 
there are reports of unresolved, steep-spectrum radio emission, consistent with jet ac¬ 
tivity, in both types of system ( Kording et al.|2008 20111. It is not clear why jets are 
not more commonly observed in these systems, although some jet launching models 
suggesting explanations exist in the literature (Soker and Lasota 2004| >. One basic 
reason may be the lack of axisymmetry in cataclysmic variables, as this is thought to 


be another major condition to launch a jet for a considerable period of time (Fendt 
and Zinnecker|l998|l. Alternatively, the observability of jet-related synchrotron emis 


sion depends strongly on the local environment of the jet, and it is possible that the 
conditions for this emission to be observed are simply not typically met in CVs. 

In this workshop proceeding we will review the many fundamental roles played 
by magnetic fields in accretion and outflows, both from a theoretical and an observa¬ 
tional point of view. The principal processes involved in jet formation can be sum¬ 
marized as follows (Blandford and Payne|1982; Pudritz and Norman|1983 Pudritz et 


al.|2007||Shang et al.|2007),' 


(1) Jets are powered by gravitational energy released through accretion and by rota¬ 
tional energy of the disk and/or the central star (or black hole). Magnetic flux is 
provided by the star-disk (or black hole-disk) system, possibly by a disk or stellar 
dynamo, or by the advection of the interstellar field. The star-disk system also 
drives an electric current. 

(2) Accreting plasma is diverted and launched as a plasma wind (from the stellar or 
disk surface) coupled to the magnetic field and accelerated magneto-centrifugally. 
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(3) Inertial forces wind up the poloidal field inducing a toroidal component. 

(4) The jet plasma is accelerated magnetically (conversion of Poynting flux). 

(5) The toroidal field tension collimates the outflow into a high-speed jet beam. 

(6) The plasma velocities subsequently exceed the speed of the magnetosonic waves. 
The super-fast magnetosonic regime is causally decoupled from the surrounding 
medium. 

(7) Where the outflow meets the ISM, a shock develops, thermalizing the jet energy. 

The overall process as outlined above is complex and takes place over a wide 
range of both length- and time-scales. Due to this complexity, the aspects of the jet 
problem have to be tackled independently. One may distinguish six principal topics 
(Fig-0- roughly corresponding to the stages in the overall picture described above. 
These include: 

(1) The accretion disk structure and its evolution, including thermal effects and the 
origin of turbulence. 

(2) The origin of the jet magnetic field, possibly through a disk dynamo, a stellar 
dynamo, or by the advection of ambient magnetic field. 

(3) The ejection of disk material into wind, thus the transition from accretion to ejec¬ 
tion. 

(4) The collimation and acceleration of ejected material into jets. 

(5) The propagation of the asymptotic jet, its stability and interaction with ambient 
medium. A related question is the feedback of jets on star or galaxy formation. 

(6) The possible presence and impact of a central spine jet, e.g. a stellar wind or 
black hole jet, in comparison to a jet originating with the disk. Under what cir¬ 
cumstances are disk jets and spine jets present or absent? 

In this article we begin by describing the theoretical support for the role of mag¬ 
netic fields in the accretion disk itself (Section |2j. We then discuss the modeling of 
magnetically driven jets (Section[3]) and their relationship to the disk (Section[4ji. Ob¬ 
servational constraints on field properties and their relation to models are discussed 
in Section[5] Our conclusions are presented in Section[6] 


2 Magnetic Fields and Accretion Disks 


We begin with the accretion disk itself, and the fundamental role played by magnetic 
fields. Throughout the history of disk theory, “shedding the angular momentum” has 
been easier said than done. The absence of an obvious mechanism to do so was con¬ 
sidered a serious objection to Laplace’s disk hypothesis for solar system formation, 


for example. By the time that the basic theory of accretion was laid out by Shakura 
|and Sunyaev| ( |1973l ) and |Lynden-Bell and Pringie] ( | 1 974| > it was accepted that Nature 
accomplished the transport of angular momentum by some means, though that means 
was unknown to astronomers. Shakura and Sunyaev introduced the “alpha viscosity” 
parameterization as a means of sidestepping the uncertainty. In this model the internal 
stress is proportional to the total pressure, Trq = aP. 

The physical viscosity in the gas within an accretion disk is known to be far too 
small to account for observed accretion rates. If the disk were turbulent, however, the 
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Fig. 1 Model sketch of MHD jet formation, indicating the six generic problems to be solved: launching, 
acceleration, and collimation; disk structure; magnetic field origin; asymptotic interaction; central spine 
jet. 


Reynolds stresses associated with that turbulence could transport angular momentum 
at the necessary rate. This seemed like a straightforward solution to the problem. Be¬ 
cause the viscosity is so small, the Reynolds number (ratio of characteristic velocity 
times a characteristic length divided by the viscosity) of the gas is correspondingly 
huge, and in terrestrial contexts high Reynolds number flows are generally turbu- 
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lent. But gas in Keplerian orbits is stable to perturbations by the Rayleigh criterion, 
which simply requires that angular momentum increase with radius, dL/dR> 0. The 
positive epicyclic frequency associated with Keplerian orbital flows is strongly sta¬ 
bilizing, and it is not energetically favorable for turbulence to develop and be sus¬ 


tained from the background angular momentum distribution (Balbus, Hawley and 
Stone|1996 Balbus and Hawley|1998| l. 


The linear stability of Keplerian flows was, of course, recognized early on as a 
difficulty for the creation of turbulence in disks. However, it was known that in the 


laboratory simple shear flows display nonlinear instability. Balbus, Hawley and Stone 


( |l996j ) argued that this nonlinear instability represented the boundary between the 
Rayleigh-stable and Rayleigh-unstable regimes; when the epicyclic frequency goes 
to zero the linear response vanishes, leaving nonlinear terms as the lowest order ef¬ 
fect. Given the importance of the question of the nonlinear stability of hydrodynamic 
Keplerian flow, several groups have carried out fluid experiments. Ji et al. ( 2006[ ) 
examined the stability of a Rayleigh-stable Couette flow and found no evidence of 
significant turbulence. Paoletti and Lathrop ( 201 1] >, on the other hand, found a break¬ 
down into turbulence in their experiment. The experiments of Schartman et al. (|2012 1 
found no transition to turbulence, and attributed the Paoletti and Lathrop ( 201 1) re- 
sults to turblence generated at the end caps of the Couette cylinder, a conclusion 
further supported by the numerical simulations of Avila (2012 1 . More recent labo¬ 
ratory experiments by jEdlund and Ji| ( j2014) ) are also consistent with stability. While 
work will no doubt continue on this important issue, at the moment the weight of 
the theoretical and experimental data lie with the nonlinear hydrodynamic stability of 
Keplerian flow. 


2.1 Basic MRI Physics 


It was, of course, recognized early on that magnetic fields could transport angular 
momentum through Maxwell stresses, —ByB^/ATZ, but unless the magnetic pressure 
were near the thermal pressure, it was expected that the effective a would be rela¬ 
tively small. As it turns out, magnetic fields render the disk linearly unstable to the 
magneto-rotational instability ( jBalbus and Hawley|l991 MRI) even when (and, in¬ 
deed, only when) the magnetic field is weak. The MRI can be visualized intuitively 
in the form of a spring connecting two orbiting masses ( jBalbus and Hawley||l998| l, 
one in a slightly lower orbit than the other. Because angular velocity decreases with 
R, Q cx R T 2 , qjg j nner mass pulls ahead of the outer mass. This causes the spring 
to stretch; the tension force pulls back on the inner mass, and pulls forward on the 
outer. The effect of this is to transfer some angular momentum from the inner to 
the outer mass, but in doing so the inner mass drops to a lower orbit, increasing its 
angular velocity, increasing the relative separation of the masses and increasing the 
spring tension force. The process runs away unless the spring is strong enough to 
force the masses to remain together, which happens only if the tension exceeds the 
orbital dynamical force. If we replace the spring with a magnetic field, the spring 
tension is replaced with the magnetic tension associated with the Maxwell stress. In 
a magnetized disk the linear stability requirement is no longer the Rayleigh criterion, 
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but is instead a requirement that angular velocity increase outward, dQ/dR > 0. A 
Keplerian disk is unstable. 

There are several remarkable properties of the MRI. One is that the question of 
stability does not depend on the magnetic field strength or orientation. The stability 
criterion. 


(k-v A ) 2 > 


d £2 2 
(ilnR 


(1) 


includes the magnetic field through the Alfven speed v \ but only in combination with 
a wavenumber vector k. Hence, for any magnetic field, however weak, one can find 
an unstable wavenumber so long as £2 decreases with R. 

The stability criterion for a Keplerian disk is 


2n 


~VA 


> V?>£2, 


( 2 ) 


^MRl 


which indicates that the unstable wavelengths, A mri, are those that are longer than the 
distance an Alfven wave travels in one orbit. The growth rate of the MRI is propor- 
tional to k ■ v A , peaking at a value of for Keplerian disks for a wavelength (Balbus 
|and Hawley|1998| l 

2 kva 


Amr/ = 


VT5 *2 


(3) 


As a practical matter, for an increasingly weak field, A mri should eventually fall to 
a scale small enough that ideal MHD no longer applies. When the MRI wavelengths 
are comparable to the resistive scale, for example, the field diffuses through the gas 
faster than the instability can grow, resulting in stabilization. At the other extreme, 
namely the strong field limit, the shortest unstable wavelength must fit into the disk. 
When A mri is comparable to the scale height of the disk, //, then ~ c s , where c s 
is the sound speed. This is a strong field indeed, and this condition is often expressed 
in terms of the plasma /3 parameter, which is the ratio of the thermal to magnetic 
pressure, as j3 < 1. Strong fields do not necessarily mean a stable, quiet disk, however. 
First, when [5 < 1 the fields are dynamically important and can transport angular 


momentum directly, possibly through a wind or a jet (see, for example, Lesur, Ferreira 
and Ogilvie ||2013). Second, the stability condition on a toroidal field requires the 


Alfven speed to be comparable to the orbital speed, rather than the sound speed; 
confining fields of this strength within a disk would seem to be problematic. 

To conclude, disks for which the gas has sufficient conductivity as to be magne¬ 
tized will be linearly unstable to the MRI. The action of the MRI is to transfer angular 
momentum outward through the disk, precisely what is needed for accretion. 


2.2 MRI-driven turbulence 

The linear analysis establishes that Keplerian orbits are unstable in the presence of a 
weak magnetic field. For any further understanding of the properties of an unstable 
disk we must turn to numerical simulations. Accretion simulations can be global or 
local. Global simulations model the whole disk, or at least a region of significant 
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radial extent. Local simulations, on the other hand, attempt to focus on the properties 
of the gas within a small region of the disk. For accretion disks, local simulations use 
an approximation known as the shearing box (Hawley, Gammie and Balbus |T995j >. 
The shearing box domain is assumed to be centered at some radial location R that has 
an orbital frequency Q (R). The size of the box is taken to be much less than R so that 
one can assume a local Cartesian geometry. The tidal gravitational and Coriolis forces 
are retained. Radial boundary conditions are established by assuming that the box is 
surrounded by identical boxes sliding past at the appropriate shear rate. Through use 
of the shearing box approximation one can study the details of MRI-driven turbulence 
on scales that are much smaller than the scale height of the disk or the radial distance 
from the central star. 

Extensive shearing box simulations have led to the following general conclusions: 


The MRI leads to MHD turbulence characterized by an anisotropic stress tensor: 
the R(j> magnetic stress component, Trq, is large, which leads to significant radial 
transport of net angular momentum (Hawley, Gammie and Balb us[1995] >. 

A net Reynolds stress is also present in the turbulence, and it too transports angu¬ 
lar momentum, but the Maxwell stress is consistently larger than the Reynolds 
stress by a factor of 3-4. If the magnetic fields are removed, the turbulence 
quickly dies out ( |Hawley, Gammie and Balhusj 1995[ ). 

Time and space variations in the turbulence can be large. The turbulence is chaotic 
(Winters, Balbus and Hawleyj |2003) l. Characterizing the stress in terms of an a 
parameter makes sense only in a space- and time-averaged sense. 

The stress is proportional to the magnetic pressure, Trq = a, nag P mag , with an 
(Xmag = 0.4-0.5. This implies that the traditional Shakura-Sunyaev value is a ~ 
1 /j3, where ji is determined by P mag , the total magnetic pressure in the turbulent 
flow. 

The magnetic energy typically saturates at a value j3 ~ 10-100. This is somewhat 
dependent on the character of the background magnetic field. Simulations with 
net vertical field tend to saturate at stronger field values. 

The local model is valid in so far as the energy released by the stresses is thermal- 
ized promptly, on an eddy turnover time-scale of order Q ~ 1 (Simon, Hawley and 
Beckwith]|2009jl. The rapid thermalization of the energy is required for a-disk 


theory to be applicable (Balbus and Papaloizo u|l999 1. 

The MRI acts like an MHD dynamo in that it sustains positive magnetic field en¬ 
ergies in the face of considerable dissipation ( [Hawley, Gammie and Balbus 1996} . 
Further, stratified shearing box simulations exhibit a behavior that can be modeled 
as an a-Q dynamo (Branden burg et al.|1995| Gressel|2010| Guan and Gammie 
201 1[>. This is, however, a small-scale dynamo; shearing box simulations have not 


provided any evidence for the existence of a dynamo capable of producing a large 
scale field. 

In simulations that include resistivity and viscosity, the amplitude of the MRI- 
driven turbulent fluctuations and the resulting angular momentum transport are 
a function of magnetic Prandtl number. Tubulence is not sustained for Pm < 1 
(Lesur and Longaretti )2607 ; Simon and Hawley 2009] ), i.e., when the resistivity 
is greater than the viscosity. 
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The local simulations tell us a considerable amount about the MRI and the re¬ 
sulting turbulence, but they cannot address questions of a global nature, such as the 
net accretion rate into the central star, losses due to winds and outflows, and the 
mechanisms by which jets might be generated. We know, however, that disks will be 
magnetized, and this implies the possibility of jet generation from the disk. In the 
next section we discuss the modeling of such jets. 


3 Magnetic Jets from Disks - Theory and Simulations 


For jet theory and simulation researchers typically consider only a subset of the issues 
listed in Section [I] and employ a variety of simplifications. An early example of a 
jet simulation that focused on the propagation of the jet and interaction with the 
surrounding medium is given by Norman et al.i( |l982j ) who carried out some of the 
first high-resolution simulations of hydrodynamic axisymmetric cylindrical jet flows. 
Some of the earliest simulations of questions related to jet launching, collimation in 
relation to disk properties were given by Uchida and Shibata] ( | 1 984[ > (see also[Shibata 
land Uchida|1985||Uchida and Shibata|1985 1. 

Since those early efforts, considerable work has been done on MHD jet modeling. 
One may distinguish i) between steady-state models and time-dependent numerical 
simulations, and also ii) between simulations considering the jet formation only from 
a fixed-in-time disk surface and simulations considering also the launching process, 
thus taking into account disk and jet evolution together. 

Clearly, many jet properties depend on the mass loading of the jet, which can only 
be inferred from a treatment of the accretion-ejection process. While numerical sim¬ 
ulations of the accretion-ejection structure potentially provide the time-evolution of 
the launching process, a number of constraints that were found, had been discovered 
already by previous steady state models (see below). 


3.1 Magneto-centrifugal Disk Winds 

Steady-state modeling of magneto-centrifugally launched disk winds have mostly 


Contopoulos 


1994a|bj). Some fully two-dimensional models have been proposed (Pelletier and Pu- 

dritz 1992 

Li 1993), including some that take into account the central stellar dipole 

(Fendt et a 


1995 

). Further, some numerical solutions have been proposed by e.g., 

Wardle anc 

Konigl (1993); Konigl et al. (2010); Salmeron et al. (20111 in weakly 


ionized accretion disks that are threaded by a large-scale magnetic field as a wind¬ 
driving accretion disk. They have studied the effects of different regimes for am- 
bipolar diffusion or Hall and Ohm diffusivity dominance in these disk. Self-similar 


steady-state models have also been applied to the jet launching domain (Wardle and 


Koni gl|1993 Ferreira and Pelletier|1993| Li|1995[ Ferreira|1997| Casse and Ferreira 


2000) connecting the collimating outflow with the accretion disk structure. In partic¬ 
ular, the fact that large scale magnetic fields need to be close to equipartition in order 
to launch jets via the Blandford-Payne mechanism is now well accepted, but was 
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first established with self-similar studies ( [Ferreira and Pelletier11993} |Li|1995[ |Fer-| 
reira |1997) . A similar comment can be made also concerning the turbulent magnetic 


diffusivity required in accretion disks. 

In addition to the steady-state approach, the magneto-centrifugal jet formation 
mechanism has been the subject of a number of time-dependent numerical studies. In 
particular, Ustyugova et al. ( 1995[ > and |Ouyed and Pudritz ( 1997[ > demonstrated the 
feasibility of the MHD self-collimation property of jets. Among these works, some 
studies have investigated artificial collimation ( Ustyugova et al.| 1999 i, a more con¬ 
sistent disk boundary condition ( Krasnopolsky et al.|1999| Anderson et al. 2005 | l, the 
effect of magnetic diffusivity on collimation ( |Fend£and Cemeljic| 2002 1, the impact 


of the disk magnetization profile on collimation ( Fendt||2006 ; 


Pudritz and Norman 


20061, or the impact of reconnection flares on the mass flux in jets from a two- 


component magnetic field consisting of a stellar dipole superposed on a disk magnetic 
field dFendt|2009j >. 

In the context of core-collapse gamma-ray bursts, Tchekhovskoy et al. (2008) 
studied the acceleration of magnetically-dominated jets confined by an external medium 
and demonstrated that jets gradually accelerate under the action of magnetic forces to 
Lorentz factors F > 1000 as they travel from the compact object to the stellar surface. 
However, Komissarov e t al,| ( [2009| ) pointed out an important problem: such confined, 
magnetized jets were too tightly collimated for their Lorentz factors to be consistent 
with observations. Does this mean that gamma-ray burst jets are not magnetically 
powered? It turns out that magnetized jets are still viable: as they exit the star and be¬ 
come deconfined, they experience an additional, substantial burst of acceleration that 
brings their properties into agreement with observations (Tchekhovskoy et al. 2010] ). 

In the context of stellar jets, Ramsey and Clarke (20111 studied the large-scale 
jet formation process spanning the whole range from the disk surface out to scales 
of more than 1000 AU. Further extensions of this model approach have included the 
implementation of radiation pressure by line forces as applied to jets in AGN (|Proga| 
et al.||2000| Proga and Kallman 2004]) or massive young stars (Vaidya et al.||201 


which may well affect acceleration and collimation of the jet material. Simulations 
of relativistic MHD jet formation were presented by Porth and Fendt (2010) find¬ 
ing collimated jets just as for the non-relativistic case (above). Applying relativistic 
polarized synchrotron radiative transport to these MHD simulation data through post¬ 
processing yields mock observations of small-scale AGN jets ( Porth et al.|2011 1 . 

The stability of the jet formation site has been studied using 3D simulations for 
the non-relativistic ( |Ouyed et al,||2003] > and the relativistic case ( |Porth|[20l3] l. Self¬ 
stabilization of the jet formation mechanism seems to be enforced by the magnetic 
“backbone” of the jet, the very inner highly magnetized axial jet region. These studies 
complement simulations investigating the stability of jet propagation on the asymp¬ 
totic scales. Here, usually a collimated jet is injected into an ambient gas distribution, 
either for the non-relativistic (Stone and Norman|1992[|Todo et al ,| 1 993[ Stone and 


Normanjj994^ Stone and Hardee|2000| l or the relativistic ( Mignone et al.|2010 Kep- 
pens et al.|2008|l case. 


In the aforementioned studies, the jet-launching accretion disk is taken into ac¬ 
count as a boundary condition, prescribing a certain mass flux or magnetic flux profile 
in the outflow. This may be a reasonable setup in order to investigate jet formation, 
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i.e. the acceleration and collimation process of a jet. However, such simulations can¬ 
not tell the efficiency of mass loading or angular momentum loss from disk to jet, or 
cannot determine which kind of disks launch jets and under which circumstances. 

It is therefore essential to extend the jet formation setup and include the launching 
process in the simulations — that is, to simulate the accretion-ejection transition. 
Clearly, this approach is computationally much more expensive. The typical time 
scales for the jet and disk region differ substantially; disk physics operates on the 
Keplerian time scale (which increases as R ' 2 ), and on the even-longer viscous and 
the diffusive time scales. The jet follows a much faster dynamical time scale. As this 
approach is also limited by spatial and time resolution, jet launching simulations to 
date have employed a rather simple disk model —namely, an a-prescription for the 
disk turbulent magnetic diffusivity and viscosity, and without considering radiative 
effects. 


Numerical simulations of the launching of MHD jets from accretion disks have 
been presented by Kudoh et al. (|1998|l; Kato et al. (2002) and Casse and Keppens 


(2002 |2004| l. These simulations treat the ejection of a collimated outflow out of an 
evolving disk through which a magnetic field is threaded. To prevent the field from 
accreting itself, the disk is resistive, allowing the field to slip through. [Zanni et~ak| 
( 2007) further developed this approach with emphasis on how resistive effects mod¬ 
ify the dynamical evolution. An additional central stellar wind was considered by 
|Meliani et al.| ( |20Q6] >. Further studies have concerned the effects of the absolute field 
strength or the field geometry, in particular investigating field strengths around and 
below equipartition (Kuwabara et al.|[2005) |Tzeferacos et al.||2009| |Murphy et al.| 


2010). These latter simulations follow several hundreds of (inner) disk orbital peri¬ 


ods, providing sufficient time evolution to also reach a (quasi) steady state for the fast 
jet flow. 

Parameter studies have been able to disentangle the effects of magnetic field 
strength (magnetization) and magnetic diffusivity (strength, scale height) on mass 
loading and jet speed (Sheikhneza mi et al.|2012) . Tzeferacos et al. (2013 1 considered 
how entropy affects the launching process, in particular how disk heating and cooling 
influence the launching process. They find that heating at the disk surface enhances 
the mass load, as predicted in the steady state modeling by Casse and Ferreira ( 2000) . 

Most jets and outflows are observed as bipolar streams. Very often, observations 
show asymmetrical jets and counterjets. For protostellar jets one exception is HH 212, 
which shows an almost perfectly symmetric bipolar structure ( Zinnecker et al.|1998) . 
For relativistic jets Doppler beaming may play a role for the observed jet asymme¬ 
tries, and, in fact, many well-behaved jets can be modeled on the assumption that 
Doppler beaming dominates the apparent asymmetry ( |Laing and Bridle|2014) . How¬ 
ever, environmental asymmetries must be also considered. It is therefore interesting 
to investigate the evolution of both hemispheres of a global jet-disk system in order to 
see whether and how a global asymmetry in the large-scale outflow can be governed 
by the disk evolution. Fendt and Sheikhnezami ( |2013) have been able to trigger jet 
asymmetries by disturbing the hemispheric symmetry of the jet-launching accretion 
disk and find mass flux or jet velocity differences between jet and counter jet of up to 
20 %. 
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von Rekowski et al.] ( [2003[ > and von Rekowski and Brandenburg ( |2004|) investi¬ 
gated the origin of the magnetic field driving the jet by including a mean-field disk 
dynamo in the star or in the disk. Asymmetric ejections of stellar wind components 
were found from offset multi-pole stellar magnetospheres (Lovelace et al. 2010} >. The 
most recent work considers the time (> 5 x 10 5 dynamical time steps) evolution on a 
large spherical grid (> 2000 inner disk radii), including the action of a mean-field disk 
dynamo that builds up the jet magnetic field. A variable dynamo action may cause 
the time-dependent ejection of jet material (Step anovs and Fendt||2014] |Stepanovs, 
|Fendt and Sheikhnezami|2014| l. 

The main jet launching processes seem to be well understood, in the sense that a 
large number of independent simulation studies give consistent results. Several issues 
are not yet resolved, however. These include: (i) the origin of jet knots, and (ii) the 
coupling of the small scale disk physics to the global jet outflow. In the end there is 
a good chance that the answers to both questions are interrelated. A full numerical 
simulation that addresses these questions would require substantial effort, including 
full 3D, high resolution, and more complete physics, in particular the treatment of 
thermal effects. First steps in this direction have already been taken. 

Concerning the origin of jet knots, it is still unclear whether knots are signatures 
of internal shocks of jet material launched episodically with different speed, exter¬ 
nal shocks of jet material with the ambient medium, or re-collimation shocks. One 
observation supporting the intrinsic origin of jet knots is that of the jet of HH212, 
which shows an almost perfectly symmetric bipolar structure with an identical knot 
separation for jet and counter jet (Zinnecker et al. |1998) . Such a structure can only 
be generated by a mechanism intrinsic to the jet source. 

The knot separation in protostellar jets typically corresponds to time scales of 
Tkin — 10 — 100 yrs. By contrast, a typical time scale for the jet launching area would 
be about 10-20 days, that is the Keplerian period close to an inner disk radius of 
0.1 AU. This time scale increases up to one year if the jet launching radius is larger, 
say up to 3-5 AU ( jFrank et al.||2014 >. However, outflows launched from such large 
radii would probably not achieve the high velocities of 300 — 500km s _1 observed 
for jets. Thus, the mechanism responsible for the jet knots must be intrinsic to the 
disk, and triggered by a physical processes on a rather long time scale. Candidates 
are disk thermal (FU Orionis) or accretion instabilities, a mid-term variation of the 
jet-launching magnetic field, or MRI-active/dead disks. 

The simulations of Stepanovs, Fendt and Sheikhnezami ( 2014j provide an exam¬ 
ple of episodic knot ejections. They apply a a-£l mean-field disk dynamo to generate 
the jet launching magnetic field. In contrary, all previous simulations of this kind start 
with a strong initial magnetic field distribution. The episodic ejections are triggered 
by switching on and off the dynamo term every 2000 dynamical time steps (Figj2|. 

The idea of a connection between the small scale disk physics and the global 
disk outflow is motivated by the fact that jet launching relies on both the existence 
of a large scale magnetic field and the existence of a (turbulent) magnetic diffusivity, 
which allows for accretion through the field lines, for angular momentum transfer, 
and, essentially, also for the mass loading onto the jet magnetic field. The disk turbu¬ 
lence is a result of the magnetorotational instability (MRI, see above) and thus of the 
small scale physics. 
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Fig. 2 Time evolution of the inflow-outflow structure for a dynamo-generated disk magnetic field. Shown 
is the density (top, colors, in logarithmic scale), and poloidal velocity (bottom, colors, linear scale), the 
poloidal magnetic flux (thin black lines), and the poloidal velocity vectors. The dynamical time step is 
about 10 5 , while the dynamo time scale (switch on/off) is 10 3 . Two ejections of material launched at 
At — 2000 are more easily seen in the velocity figure (below), than in the top figure, which shows a 
slightly later time when the preceding knot has left the grid already. Length units is the inner disk radius 
of ~ 0.1 AU. Figure taken from |Stepanovs and Fendt|j2014) . 
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First results combining a mixed numerical-analytical treatment indicate that large 
scale MRI modes may produce magnetically driven outflows, however, these flows 
seem to be 3D-unstable (Kuwabara et al. 2005| , and thus may be a transient ef¬ 
fect. On the other hand, local numerical simulations suggest that magnetocentrifugal 



have found that the MRI near equipartition does not lead to turbulence, but can be 
responsible for jet launching. 

The overall goal is still to answer the question: what kind of disks launch jets and 
what kind of disks do not? It is clear from statistical arguments for both protostars and 
extragalactic sources that jets are a relatively short-lived phenomenon. The typical jet 
propagation time scale for protostellar jets is about Td yn = Tjet/Fjet — 10 4 yrs, while 
protostellar disks live some 10 6 yrs. A similar argument can be raised for AGN jets 
based on the fact that there are far more radio-quiet AGN (quasars) than radio-loud 
AGN (jet sources), although at least some types of radio-loud AGN appear to be 
almost continuously in an ‘on’ state. So far, modeling has not provided an answer to 
this question. 


4 Black Hole Accretion, Fields and Jets 


4.1 Simulating the black hole jet 


As we have seen, a plethora of models and simulations have demonstrated the ability 
of magnetic fields to launch, accelerate and collimate jet outflows. But simulations 
designed to study jets have, as discussed in the previous section, typically used one or 
another specific set of initial or boundary conditions that pre-suppose a disk stmcture, 
a mass-energy flux, a magnetic field configuration, or some combination of these. An 
alternative approach is to focus not on the jet, but on on the accretion disk itself and 
see when and under what circumstances a jet might develop from the flow. 

Most global simulations begin with an orbiting axisymmetric hydrostatic equilib¬ 
rium (torus) a few tens of gravitational radii from the black hole. The initial magnetic 
field is entirely contained within the matter, so that there is no net magnetic flux 
and no magnetic field on either the outer boundary or the event horizon. A favorite 
configuration consists of large concentric dipolar loops (Hawley||2000[ McKinney 


and^ummie|2()()4j JDeJVilhers, Hawley, Krolik and Hirose|2005 [Hawley and Kro 


lik|2006 McKinney ||2006 1 . The development of the MRI subsequently leads to the 

establishment of outflows from the disk. The wind generation seems to be largely an 
outcome of thermodynamics. In simulations that do not include cooling, the combi¬ 
nation of the heat released from accretion and the buildup of magnetic pressure lead 
to pressure-driven outflows. In the absence of a large scale vertical field through the 
disk, however, this wind is neither collimated nor unbound. 

Global simulations have nevertheless produced jets that are consistent with the 
Blandford-Znajek mechanism, through the creation of a substantial axial field that 
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threads the black hole. For example, in a model that begins with a toms containing 
dipolar field loops, differential rotation rapidly generates toroidal field. This field in¬ 
creases the total pressure within the toms, driving the its inner edge inward. As the 
field is stretched radially, the Maxwell stress transports angular momentum outward, 
further enhancing the inflow. Subsequent evolution carries the inner edge of the disk, 
and the magnetic field within, down to the black hole. The field expands rapidly into 
the nearly empty funnel region along the black hole axis. Gas drains off the field lines 
and into the hole leaving behind a low axial field. Frame dragging by the rotating 
hole powers an outgoing Poynting flux. The power of the Poynting flux jet is deter¬ 


mined by the spin of the black hole and the strength of the magnetic field (Blandford 


and Znajek [l977) , as we discuss below. The simulated jet’s power matches well with 


the predictions of the Blandford-Znajek model (McKinney and Gammi ej2004[|McK-| 
inney|2005j Tchekhovskoy et al.|201 1[ >. 


The matter content within the Blandford-Znajek jet is very small; the angular 
momentum in the disk material is too great to enter the axial funnel. In the simulations 
the jet density is typically set by the value assigned to the numerical vacuum. How 
these jets become mass-loaded in Nature is still something of an open issue. Since 
the boost factor of the jet is determined by the mass loading, F is typically large, e.g., 
F ~ 10 ( |McKinney|2005j l. No particular weight can be given to any value of F found 
in simulations, although it is safe to say that the feasibility of values as large as those 
observed (Section]?]) has been demonstrated. 

Although there is no appreciable matter within the funnel, some of the field lines 
that lie just outside the funnel pass through the accretion flow, and those field lines can 
accelerate matter outward in a collimating flow. This “funnel wall jet” was observed 
in early pseudo-Newtonian simulations (Hawley |2002)>, and is a regular feature of 
fully relativistic jet simulations, e.g., Hawley and Krolik ( 2Q06| >. 

Because the field lines in the jet rotate, the jet carries angular momentum as well 
as energy away from the black hole. The angular momentum flux can be large, com¬ 


parable to that brought down by accretion ([Gammie, Shapiro and McKinney |2004 


Hawley and Krolik|2006[ Beckwith, Hawley and Krolik 2008). One potentially im¬ 


portant consequence of this was noted by Gammie, Shapiro and McKinney| ( |2004| ) 
who pointed out that electromagnetic angular momentum flux in the jet rises so 


steeply with black hole spin that it may limit a/M to ~ 0.93 (see also 

McKinney 

and Gammie 2004, Hawley and Krolik 2006 Beckwith, Hawley and Krolik 2008) or 

even down to < 0.1, as we discuss below ( 

Tchekhovskoy et al. 2012a 

Tchekhovskoy 


|2015| ). The possibility that magnetic fields could play such an important role in black 
hole physics is certainly in keeping with the spirit of this workshop! 

For a magnetic field strength B at the event horizon, jet power is (roughly) given 
by the magnetic energy density, B 2 , times jet cross-section at its base, /" (where 
r g = GM/c 2 is the black hole gravitational radius and M its mass), times the speed 
v ~ c at which the field moves out. 


Fje t ~ ( a/M) 2 B 2 r g c ~ (, a 2 c/M 2 r 2 g )<P 2 = k<P 2 , (4) 

where the effect of black hole rotation is in the ( a/M ) 2 pre-factor, and the magnetic 
field is expressed in terms of the flux <J> ~ Br 2 . Since in the course of an observation 
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Fig. 3 Panels (a) and (b) show a meridional slice of density with color (red shows high, blue low den¬ 
sity, see color bar) and magnetic field lines with black lines for the initial conditions of the simulations 
of |Tchekhovskoy et al.|(20TT| > The initial magnetic field is fully contained within the torus, and have a 
sufficient amount of magnetic flux to saturate the black hole with magnetic flux and to lead to the develop¬ 
ment of a magnetically-arrested disk (MAD). Panel (c) contrasts this initial condition with a typical initial 
condition in earlier simulations where the torus extent is smaller, and the magnetic flux is reduced. This 
is labelled “SANE” ICs, for “standard and normal evolution” {Narayan et al.|2012). Figure adapted from 
|Tchekhovskoy|{20l5) . 


M and a usually do not change appreciably, the factor k = a 2 c/M 2 r 2 is a constant, 
and all changes in jet power are due to changes in the magnetic flux, <f>. 

What is the possible range of <J> and the corresponding range of ^jet? The lower 
limit is clear: zero magnetic flux leads to zero power. But what sets the maximum 


value of <P ? In the simulations of Beckwith, Hawley and Krolik (20091 the black 
hole field strength is comparable to the (mostly gas) pressure in the inner disk, /3 ~ 1, 
suggesting that, in general, equipartition holds. This further suggests that even higher 
black hole magnetic field strengths are possible for more energetic accretion flows. 
Tchekhovskoy, Narayan and McKinney (2011|) carried out numerical simulations of 


black hole disk-jet systems designed specifically to address the question of maximum 
field strength. The simulation begins with an initial thick, hot, large scale torus (as 
shown in Fig. [3ji,b). Because of the large torus size, it contained a particularly large 
amount of magnetic flux, much larger than in previous work. Fig. |3ji,b shows the 
initial condition, compared with initial conditions typically used in previous simula¬ 
tions (Fig. [3p). Using these initial conditions Tchekhovskoy et al. ( 2011 1 ) found that 
magnetic flux on the black hole saturated at an equipartition point, one where the 
magnetic pressure was sufficient to halt accretion against the inward pull of gravity. 
In other words, the dynamically-important magnetic flux obstructed the accretion and 
led to what is dubbed a magnetically arrested disk, or a “MAD”; in such a condition 
no further flux could be carried down to the hole. In contrast, models where magnetic 
field falls short of saturating the black hole were referred to as “SANE” accretion 
models, for “standard and normal evolution” by Narayan et al. (2012 1 . In the MAD 


simulation <f> and Pj et are greatly increased, as we discuss below (Tchekhovskoy et 
|al.|2011||20l2al|McKinney et al.|2012||Tchekhovskoy|20151 l. 
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How do such strong magnetic fields manage to stay in force-balance with the 


disk? As in previous disk simulations, e.g., Beckwith, Hawley and Krolik(2009 1 , the 
black hole field strength is comparable to the total pressure in the disk. In MADs, 
magnetic and gas pressures in the disk are comparable, J3 ~ 1, but the much stronger 
black hole magnetic field in the MAD case compresses the inner disk substantially, 
so the inner regions of MAD disks are actually over-pressured compared to previous 
simulations by a factor of ~ 5 —10. 

Figure [4] shows the structure of the simulated disk-jet system. An accretion disk 
around the central black hole leads to a pair of jets that extend out to much larger 
distances than the black hole horizon radius and are collimated into small opening 
angles by outflows from the disk (not shown). Jet magnetic field is predominantly 
toroidal, which reflects the fact that the jets are produced by the rotation of the black 
hole, which twists the magnetic field lines into tightly-wound helices. 

Even though the accretion disk and the jets turn out to be highly time-variable, 
in a time-average sense the structure of the flow is remarkably simple. The poloidal 
magnetic field lines are shown with black solid lines in Fig.[4jb). The group of field 
lines highlighted in green connects to the black hole and makes up the twin polar 
jets. As discussed above, these field lines have little to no gas attached to them; disk 
gas cannot cross magnetic field lines and has too much angular momentum to enter 
the funnel. The large Poynting flux and very small inertia yields highly relativistic 
velocities. The field lines highlighted in blue connect to the disk body and make up 
the magnetic field bundle that produces the slow, heavy disk wind, whose power is 
much smaller than P }et for rapidly spinning black holes. 

The jet power can be characterized by the outflow energy efficiency rj as 


V = 


^outflow 

Me 2 


(5) 


For the MAD simulation, the accretion rate M fluctuates considerably, but does main¬ 
tain a well-defined average value. In the MAD flow, for a rapidly spinning black hole, 
a/M = 0.99, the total outflow power, P 0 utflow> which is the sum of the black hole 
powered jet, /jet and the disk powered wind, exceeds the accretion power, 77 > 1, 
i.e., net energy is extracted from the black hole (Tchekhovskoy et al'| 20l 1 [ ). This is 
a promising result since high values of jet power are consistent with or are preferred 
by observations (RawlingsandSaunders|1991| Ghisellini et al. 201 0[ Fernandes et 


al.|201 1| |McNamara et al.|201 1[ Nemmen and Tchekhovskoy|2014[ Ghisellini et al. 

20141). 


The MAD flow hearkens back to the ion-torus model of |Rees et al.| ( |1982| where 
the role of the thick torus is to confine and anchor the fields that extract rotational 
energy from the hole, and accretion per se is near zero. This model was put originally 
forward as an explanation for systems where high jet power is seen, but very little 
luminosity emerges from the central engine. In the MAD flow, accretion continues at 
a non-zero rate regulated by the magnetic interchange instability in the disk and the 
strength of the central magnetic field. 

Figure|5 compares an analytic approximation for /j et in MAD simulations to some 
earlier simulations (McKinney 2005] [Hawley and Kro lik 2006). Jet efficiency in the 
MAD case is a 2 , suggesting the ratio of magnetic flux squared to mass accretion 
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Fig. 4 [Panel (a)]: A 3D rendering of a MAD disk-jet simulation. Dynamically-important magnetic fields 
are twisted by the rotation of a black hole (too small to be seen in the image) at the center of an accretion 
disk. The azimuthal magnetic field component clearly dominates the jet structure. Density is shown with 
color: disk body is shown with yellow and jets with cyan-blue color; jet magnetic field lines are cyan 
bands. The image size is approximately 300r^ x 800r 5 . [Panel (b)]: Vertical slice through a MAD disk- 
jet simulation averaged in time and azimuth. Ordered, dynamically-important magnetic fields remove the 
angular momentum from the accreting gas even as they obstruct its infall onto a rapidly spinning black 
hole ( a/M — 0.99). Gray filled circle shows the black hole, black solid lines show poloidal magnetic field 
lines, and gray dashed lines indicate density scale height of the accretion flow, which becomes strongly 
compressed vertically by black hole magnetic field near the event horizon. The symmetry of the time- 
average magnetic flux surfaces is broken, due to long-term fluctuations in the accretion flow. Streamlines 
of velocity are depicted both as thin red lines and with colored “iron filings”, which are better at indicating 
the fine details of the flow structure. The flow pattern is a standard hourglass shape: equatorial disk inflow at 
low latitudes, which turns around and forms a disk wind outflow (labeled as “disk body” and “disk wind”, 
respectively, and highlighted in blue), and twin polar jets at high latitudes (labeled as “jet” and highlighted 
in green). Figure taken from |Tchekhovskoy|{2015) . 


rate is roughly constant ( Tchekhovskoy et al.|2012aj Tchekhovskoy|20f5] l. If MADs 
achieve the maximum possible amount of magnetic flux threading the black hole then 
they have achieved the maximum possible jet power for a given mass accretion rate 


(see Fig. pi. The other models show a rapid rise as a/M approaches 1. Hawley and 


KrolikI ( |2006| > attribute this to greater local field amplification for the extreme Kerr 
holes; the field has “room to grow,” and is a factor of a few smaller than in the MAD 
simulations. 

Efficiencies T] > 1 can lead to spindown of the black hole. For example, in an 
AGN accreting at 10% Eddington luminosity, the central black hole is sp un down 
to near-zero spin, a/M <0.1, in x ~ 3 x 10 8 years (Tchekhovskoy et al. 2012a 


|Tchekhovs koy 2015). This value is interesting astrophysically because it is compa 
rable to a characteristic quasar lifetime (Section [5]). Over this time period, jets can 
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Fig. 5 Comparison of jet energy efficiency obtained in MAD simulations, Pj a /Mc 2 , which is shown 
with green solid line, with previously reported approximations of simulated jet power: Hawley and Krolik 
{2006} , which is shown with black dash-dotted line (HK06), and McKinney (2005), which is shown with 
brown dashed line (M05), plotted over the range 0 <a/M < 0.99. MADs can produce much more powerful 
jets because they have the maximum possible amount of magnetic flux threading the black hole. 


extract a substantial fraction of the central black hole spin energy and deposit it into 
the ambient medium. The central galaxy in the cluster MS0735.6+7421 may be one 
such example (|McNamara et al. 12009). 


4.2 Origin of the jet field 

Jets appear to require some large-scale magnetic field, in particular a net vertical field 
connected either to the disk itself or the central star or black hole. This raises the 
question of the origin of that field. In the jet-forming global simulations, the initial 
condition consists of dipole loops embedded within initial orbiting gas, typically a 
torus of some thickness. In essence, the subsequent evolution (which is a consequence 
of the initial conditions) transfers the dipole field down to the black hole. Other initial 
field configurations, such as quadrupolar, toroidal, smaller scale loops, have been 
found to be far less effective in creating jets ( jBeckwith, Hawley and Krolik|2008j >. 

Simulations that produce jets rely on favorable initial conditions. What happens 
in Nature? At least two possibilities exist. One is that the required field is gener¬ 
ated in situ by the disk itself through a dynamo process, as a natural consequence of 
accretion. Some example dynamo models were discussed above (e.g.. Fig. [2]). The 
other possibility is that the required held is carried into the central disk/star, from 
the companion star in a binary system and from the interstellar medium in an AGN. 
The advantage of advection is that a relatively weak net held at large radius could be 
signihcantly amplihed simply by the geometric factor in going from large to small 
radius. For example, magnetic helds at the edge of the sphere of influence of a su- 
permassive black hole, i.e., r ~ 100 pc, are plausibly B ~ jj. G. Suppose these helds 
maintain their coherence over roughly a similar length scale. The magnetic flux con¬ 
tained in patch of held this size, <f> p;uc h ~ 10 35 G cm 2 , is much larger than the flux 
necessary to saturate a black hole with magnetic held, <J>mad ~ 10 33 ' 5 G cm 2 (for a 
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black hole of mass M = 10 9 M Q accreting at 10% of Eddington luminosity Narayan 
et al.|2003l Tchekhovskoy|2015] l. 

Under what general circumstances can a net field be advected in by the accretion 
flow? The question is traditionally framed in the terms used by van-Ballegooijen 
( 1989) 1, namely that accretion rate is set by the viscosity v and the magnetic diffusion 
by the resistivity T] . Whether or not a field can be advected inward depends on the 
magnetic Prandtl number P m v/rj. If P m exceeds 1 then field advection will not 
occur (see also Lubow et al. 1994] ). In accretion disks the viscosity and resistivity are 
not the physical values, but instead are effective values resulting from the turbulence 
( |Guan and Gammie|2009||Lesur and Longaretti|2009[|Fromang and Stone|2009| l. 

In global simulations, flux in the initial torus (e.g. Fig. [5]) finds its way onto the 
black hole, despite starting from a large radius. This suggests that geometrically- 
thick disks are capable of transporting large-scale magnetic flux inward, at least from 
the distances of r ~ 1 ()()/-„. In standard a-disk theory, infall rates are proportional to 
(, h/r ) 2 . If h/r is close to unity the infall time becomes relatively rapid, and the inflow 
rate can exceed the diffusion rate. Thick, radiatively inefficient flows, then, seem to 
have the best chance of transporting flux inward. 


The question of whether a turbulent, radiatively efficient disk could transport flux 
inward remains an interesting one. Beckwith, Hawley and Krolik ( 20Q9) > studied the 
question through a global simulation of an initial orbiting torus with an h/r ~ 0.1 em¬ 
bedded in a vertical magnetic field. They found that there was no net flux transported 
through the disk itself; the turbulence transported angular momentum and mass, but 
not net flux. Net flux was nevertheless brought down to the black hole by a coronal 
inflow surrounding the accretion disk. Initially the field lies along cylinders of con¬ 
stant rotation, but at high altitude B 2 /p is much larger than deep inside the disk, and 
the field is subject to a finite-amplitude version of the MRI. The sign relation between 
the radial field component and the toroidal field component is the usual one, so the 
Maxwell stress has the right sign to transport angular momentum outward. The angu¬ 
lar momentum flux due to the Maxwell stress is large compared to the fluid angular 
momentum density. Fluid elements at high altitude, both above and below the equa¬ 
tor, are therefore driven inward quickly, carrying half-loops of field with them. When 
one loop approaches the central black hole from above and outward while another 
approaches the disk from below and outward, their local field directions are opposed, 
and they can reconnect. Reconnection changes the field topology, creating a closed 
loop at larger radius and an open field line at small radius. The open field line, the 
one carrying the net flux, is located where the reconnection occurred, which is at a 
radius considerably smaller than the initial radius of the field line. Thus, reconnection 
causes the flux to move inward in large jumps. All of this takes place faster than the 
mean mass inflow rate within the disk at the vicinity of the flux line’s initial position. 


In a similar vein, [Suzuki and Inutsuka[|2014| carried out a series of simulations 
with a relatively thick disk (sound speed ~ 0.1 of the orbital speed) with an initial 
weak magnetic field. In these simulations net vertical flux was carried inward by a 
rapdily infalling layer near the disks surface. All these simluations transport net flux 
through a mechanism that is governed by large-scale torques and resulting rapid in¬ 
flows, rather than by small-scale turbulence where the turbulent velocities exceed the 
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mean inward drift velocity associated with accretion. It appears plausible, therefore, 
that net inward flux transport can occur for certain types of accretion flows. 

Testing jet-generation models, including quantitative predictions for the magnetic 
flux close to the black hole, fundamentally depends on observations that can give 
insights into the field properties both at the jet base and further out. We discuss those 
in the following section. 


5 Observational constraints on jet magnetic fields 

5.1 Introduction: basic jet properties 


In this section we ask what the observational constraints are on the field strengths and 
configurations in jets, and whether these observational results can be connected to jet¬ 


launching models. (For a more detailed review of these issues see Pudritz, Hardcastle 


and Gabuzda (2012); parts of that work are summarized in this section.) We focus on 


relativistic AGN jets, since these are by far the best-studied of the systems producing 
synchrotron radiation, and begin by setting out some basic properties of the jets in 
radio-loud AGN. 

AGN jets are fast outflows carrying kinetic^ powers that are estimated to be be¬ 
tween ~ 10 42 and ~ 10 47 erg s “ 1 (e.g., Rawlings and Saunders 1991). Bulk speeds 
are known to be highly relativistic on small scales (F ~ 10-30, and perhaps higher 
in some objects: Lister et al. 2009a ; on kpc scales some jets have decelerated to sub- 
relativistic speeds ( Laing and Bridlej2002 1 while in more powerful systems jets are 
at least mildly relativistic out to hundred-kpc scales ( Mullin and Hardcastle|2009] ) 
and may retain a highly relativistic spine with F ~ 10 out to those scales, though 
direct observational evidence for this remains debatable (e.g., Tavecchio et al.|2000| 


|Hardcastle|2006| [Marshall et al.|20TT]|Konar and Hardcastle|2013[l. 

Jet composition remains an open question. Observations of synchrotron radiation 
require relativistic leptons and magnetic fields to be present, but direct observational 
constraints on the presence of a hadronic component are hard to obtain. Some indirect 
constraints suggest that the bulk of the required energy for jets can be carried by the 
radiating leptons without any need for additional components on scales from pc to 
kpc (e.g., Wardle et al.|1998 Croston et al. 2005 [ Wykes et al.|2013), It is certain that 


all jets must entrain material from stellar winds in the host galaxy (B icknell|1994| l so 
their composition may change with distance along the jet, and there is some evidence 


from pressure balance arguments that this is indeed the case (Croston and Hardcastle 
[20141 ). 


One important observational constraint on jet launching mechanisms comes from 
the fact that they appear to be able to originate in AGN with a wide range of different 
properties. The vast majority of radio-loud AGN are so-called ‘low-excitation radio 
galaxies’ which show no evidence at any wavelength for a standard radiatively effi¬ 
cient accretion disk ( [Hardcastle, Evans and Croston|2009~[|Best and Heckman[2012[ ). 
However, the existence of powerful radio-loud quasars with L/L\ m ~ 1 in the optical 
shows that radiatively inefficient accretion is not directly connected to radio jet ac- 
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tivity. In fact it appears that radio-loud AGN have a very wide range of L/Lem, even 
when jet power is taken into account ( Mingo et al.|[2014) >. 

Finally, it is important to note that observations of synchrotron radiation require a 
particle acceleration mechanism. There is indirect or direct evidence for local particle 
acceleration on all scales of radio-loud AGN up to scales of tens of kpc. Shocks are 
almost certainly implicated — some of the best-understood sites of particle accelera¬ 
tion, the hotspots of powerful double radio galaxies, are clearly physically associated 
with the jet-termination shock and show properties consistent with simple first-order 
Fermi acceleration fMeisenheimer, Yates and Roser[|1997| . However, other mecha¬ 
nisms may well be necessary to explain diffuse particle acceleration distributed on 
scales of many kpc (e.g., Hardcastle, Croston and Kraft||2007] >. When interpreting 
observations, it is necessary to bear in mind that the particles being observed at a 
given location may owe their energy spectrum to acceleration in some region that is 
spatially quite distinct. 

Because jets can be persistent (ages of the order 10 8 years have been estimated 
for some sources based on spectral or dynamical age estimates) they produce large- 
scale, long-lasting structures consisting of material that has passed up the relativistic 
jet (Scheuer 1974| , generally called lobes or plumes. The physical conditions in these 
structures are to some extent determined by the jet and so they give important clues 
to aspects of the jets that are hard to investigate directly. Because the observational 
techniques used are rather different, it is usual to discuss the large-scale lobes, the 
kpc-scale jets and the pc-scale jets separately, a convention that we follow in the 
remainder of this section. 


5.2 Observational tools 

The combination of relativistic leptons and magnetic field gives rise to synchrotron 
emission, while relativistic leptons and a photon field give us inverse-Compton emis¬ 
sion. These are the basic tools available to estimate the magnetic field properties in 
jets and their products. 

Optically thin synchrotron radiation cannot on its own give a measurement of 
magnetic field strength. For simplicity, let us represent the electron energy spectrum 
as a power law with index p, so that N(E ) = N fi E p . Then the synchrotron emissivity 
in the optically thin regime is given by 

(p+1) 

J(v) = CN 0 v —2 8 * ( 6 ) 

where C is a constant depending only weakly on p. Equation[6]implies that the same 
observed emissivity can be produced by any combination of the number density of 
electrons (scaling as No) and the strength of the field B. However, total intensity emis¬ 
sion does give important information about the magnitude of any spatial variation in 
field strength (combined with electron density, energy spectrum etc) given the strong 
dependence on B of emissivity (eq. [6]». The direction of observed polarization in opti¬ 
cally thin regions gives us an (emission-weighted, line-of-sight integrated, projected) 
estimate of the local direction of magnetic field, which is our best probe of the vector 
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properties of the field, while fractional polarization tells us about its ordering: for 
a uniform field the fractional polarization for power-law electron energy index p is 

n = (p+i)/(p + 3 ). 

Faraday rotation can complicate the interpretation of polarization. For the simple 
case of an external Faraday screen, the rotation angle (j) is given by 



(7) 


where K is a constant, n t h is the number density of thermal electrons, and the integral 
is along the line of sight to the source. If the value of (j) varies within the resolution 
element, depolarization rather than simple rotation with (j) A 2 will be observed. In 
some cases, therefore polarized intensity, and in particular the dependence of frac¬ 
tional polarization or position angle on frequency, actually tells us more about the 
foreground magnetoionic medium (in the host or the Milky Way) than it does about 
the source itself. The effects of internal and external Faraday-active media on po¬ 
larization are well understood in some simple analytic cases (jBum|1966| Jones and 
|Q’ Dell| 1977 [ |Cioffi and Jones|1980) . 


Optically thick synchrotron radiation is an essential tool for field strength mea¬ 
surements on the pc scale; the turnover frequency for synchrotron self-absorption is 
sensitively dependent on the magnetic field strength, with the absorption coefficient 
going as #(p+ 2 )/ 2 . However, as the turnover also depends on the normalization of 
the electron energy spectrum, this method depends on good information about the 
structure of the components where the turnover is observed. 

The combination of inverse-Compton and synchrotron emission gives an excel¬ 
lent constraint on the mean field strength, if both processes can be measured from 
the same region; inverse-Compton emissivity depends only on the number density of 
electrons (for a known photon field, often the CMB) so that the field strength may be 
directly estimated from the synchrotron emissivity. Again, it is necessary to know the 
geometry of the emitting region accurately. 

Finally, if all else fails, it is traditional to resort to the assumptions of equiparti- 
tion. 


Vn= ! EN{E)iE 


( 8 ) 


or minimum energy ( Burbidge|1956 1 


w “' = £ + /™ (£)d£ : “sf ' 0 <9) 

which, putting in B-dependent estimates of No from eq. [6] give fairly similar mag¬ 
netic field values for a given volume emissivity. Until recently there has been little 
observational justification for the equipartition assumption (though see below) and it 
should still be applied with caution. 
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5.3 Observations: large-scale components 


In the large-scale lobes and plumes of radio galaxies the gold standard for field 
strength measurements is the inverse-Compton technique. Large numbers of power¬ 
ful (FRII: Fanaroff and Riley 1 91 A) radio ga laxies now have global lobe field strength 
measurements (Kataoka and Stawarz 2005; Croston et al. 2005) using X-ray inverse- 
Compton measurements made with XMM and Chandra. These imply field strengths 
close to, but somewhat below, the equipartition values: a typical field strength in a 
lobe is of the order 10 jJ. G. Fewer such measurements exist for the lower-power FRI 
radio galaxies, in part because of their typically brighter thermal environments but 
also because some part of their internal pressure is probably provided by thermal 
material (Hardcastle and Croston 2010| . However, a magnetic field strength around 
1 jJ. G, again close to the equipartition value, has been estimated for the nearby FRI 
radio galaxy Centaurus A based on Fermi y-ray detections of its large-scale lobes 
dAbdo et al.|20T0] >. 

Inverse-Compton observations also provide important evidence for the spatial in- 
termittency of the magnetic field: the inverse-Compton surface brightness (which, for 
scattering of the CMB, depends only on the electron energy density and the line-of- 
sight depth) is seen to be much more uniform than the synchrotron surface brightness 
(Hardcastle and Croston f2005 1. This is as expected if the lobes are turbulent with 
an energetically sub-dominant magnetic field, and similar effects are seen in numer¬ 
ical MHD modeling of lobes (Tregillis, Jones and Ryu||2004| Hardcastle and Krause 
f20l4j ). 

Turning to field direction, the apparent field directions inferred from polarization 
measurements in the lobes are almost always perpendicular to the jet far away from 
the jet termination, with mean fractional polarization in the range 20^-0% at high 
frequencies. This is consistent with the idea that a toroidal component of the field 
dominates on large scales, as would be expected if an initially disordered field ex¬ 
pands into the lobes, but compression of a disordered field can give similar results 
( |Laing|1980| . 


5.4 Observations: kpc-scale jets and termination features 


Inverse-Compton measurements of field strengths in jets are difficult because of the 
existence of in situ particle acceleration, which can give rise to high-energy (e.g. 
X-ray) synchrotron emission. Arguments for high bulk speeds in kpc-scale jets, dis¬ 
cussed above, assume that the X-ray emission from these jets is synchrotron, but 
there is a degeneracy between the magnetic field strength and the (a priori unknown) 
jet speed and angle to the line of sight, so that these observations are not very useful 
for measurements. In low-power jets, where the dynamics are better understood, TeV 
y-ray observations place some constraints on magnetic field strength (Hardcastle and 
Croston 2011] ), implying B > 20 /iG to avoid over-producing the observed TeV emis¬ 
sion. These models depend on a good knowledge of the photon field, in this case the 
starlight at the centre of the elliptical host galaxy. 
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The magnetic field direction is much better understood. In powerful jets the field 
direction is almost universally along the jet (although it is not clear whether the syn¬ 
chrotron emission from the jet is a reliable tracer of the jet itself — it may simply 
show a boundary layer). In low-power jets, observations of polarization are essen¬ 
tial to models that break the degeneracy between jet speed and angle to the line of 
sight ( Laing and Bridle||2002| ) and when this is done it is found that the field tends 
to evolve from a predominantly longitudinal configuration on the smallest scale to 
a predominantly toroidal one after jet deceleration: no significant radial component 
is required by the models (Laing and Bridle 2Q14| >. Helical field models on these 
scales for these well-studied low-power jets are ruled out by the lack of transverse 
polarization asymmetry. 

Powerful jets terminate in compact (1-10 kpc scale), bright features known as 
hotspots, which are thought to trace the jet termination shock. As the downstream 
material should be subsonic, and so at most mildly relativistic, hotspots are easier 
to understand than the jets themselves. Their brightness means that they are good 
sources of inverse-Compton emission via the synchrotron self-Compton process, pro¬ 
ducing X-ray emission (Harris, Carilli and Perle yfT994} . Inverse-Compton modeling 
implies field strengths of the order 100 jJ .G in hotspots of powerful sources (Hard 


[castle, Harris, Worrall and Birkinshaw|2004| Kataoka and Stawarz|2005 ). Magnetic 


field structures in hotspots, as revealed by high-resolution polarization imaging, are 
complex, but there is a general tendency for the projected field direction to be per¬ 
pendicular to the jet (Leahy et al.|1997|. 


5.5 Observations: pc-scale jets 


Inverse-Compton methods are very uncertain on small scales—high-energy emission 
from the jet is often hard to distinguish from emission from other components of 
AGN. In objects where this is not the case, such as blazars, the jet geometry, bulk 
speed and angle to the line of sight are likely to be poorly known. Most methods 
for estimating the magnetic field strength on parsec scales rely on inferences from 
very long baseline interferometry (VLBI) in the radio, and hence on synchrotron 
emission. Typically, the parsec-scale radio morphology of a bright AGN manifests 
a one-sided jet structure due to Doppler boosting, which enhances the emission of 
the approaching jet. As noted above, Lorentz factors on these scales can be high, 
r ~ 10-30: constraints come from observations of apparent motions (Cohen et al. 


2007 


20091 


Lister et al.|201 3|> and from radio variability (|Jorstad et al.|2005||Hovatta et al. 


S_avolainen et a l,12010) 1. 


The apparent base of the jet, which is often the brightest and least-resolved part, 
is commonly called the core. As noted above, synchrotron self-absorption, which de¬ 
pends on both the particle number density N e and the magnetic field B. will affect 
the appearance of the jet; the core is thought to represent the jet region where the 
optical depth is equal to unity. As the number density and magnetic field depend on 
distance from the jet base, the position of the apparent core is expected to depend on 
observing frequency. This is the basis of the core-shift method for estimating mag¬ 
netic field strengths (Lobanovj 1998] [O’ Sullivan and Gabuzda||2009b] l, evaluating the 
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relation between core flux at given frequency and mass of black hole (Heinz and 
Sunyaev|2003|l, jet composition ( Hirotani|2005| l, and flow magnetization (Nokhrina 


et al.||2015|l. This method necessarily assumes that the jet is well described by the 


simple, homogeneous conical jet models of Blandford and Konigl ( 1979j>. The lat¬ 
ter assumption is supported by multi-frequency observations by (Sokolovsky et al. 

f20TT ). 

As an alternative to the core-shift method, individual components of the jet can 


be fitted with self-absorption models (e.g. ISavolainen, Wiik, Valtaoja and Tornikoski 


20081; this method is less model-dependent but has larger random errors as a re¬ 


sult of the uncertainties on component size, geometry and Doppler factor. Magnetic 
field strengths derived by either method are of the order 1 G on pc scales. Recently, 
|Zdziarski et al.| ( ]2014) > have proposed a “hybrid” method of doing core-shift analysis 
which does not rely on the equipartition assumption, again obtaining similar magnetic 
field strength values. 

On pc scales, the inferred jet magnetic field direction is normally either paral¬ 
lel or perpendicular to the jet direction (e.g. Lister and Homanj |2005| >, with jets in 
more powerful objects tending to have parallel fields. Misaligned polarization that 
is definitively not the result of Faraday rotation is rare. The interpretation of these 
field structures varies: they are often attributed to internal jet dynamics, e.g. shocks 
in the case of transverse apparent field direction, but could equally well be the result 
of an underlying helical magnetic field geometry with varying pitch angle (Lyutikov, 
Pariev and Gahuzda|2005). (As noted above, helical fields are expected in most jet 


generation models due to the toroidal component introduced by field winding: see 
Sections 1|4.2 ) As pointed out by Blandford ( |1993j >, the helical-field model would be 
supported by observations of transverse gradients of Faraday rotation across the jet, 
provided that there are thermal electrons in the jet or in its immediate environment 
to provide the required magnetoionic medium (eq. [7]). Detection of such gradients 
is observationally difficult, since the jets are poorly resolved, but they have been ob¬ 
served in a number of objects (e.g., Asada et al.|2002[[0’Sullivan and Gabuzda 2009a 


Croke, O’Sullivan and Gabuzda |2010 . While their interpretation is still controversial 
( [Tayl or and Zav ala|2010 1 these gradients provide at least some direct evidence for he¬ 
lical field structures in the inner parts of the jets. 


5.6 Core-shifts and jet magnetization 

Two important MHD parameters for describing relativistic flows are the Michel mag¬ 
netization parameter Om and the multiplicity parameter A. The first one tells us how 
strongly the flow is magnetized at its origin, and it determines the the maximum 
possible bulk Lorentz factor of the flow. The second one is the dimensionless multi¬ 
plicity parameter A = /V e //Vc,j, which is defined as the ratio of the number density N e 
to the Goldreich-Julian (GJ) number density /Vgj = £213/2k ce — the minimum num¬ 
ber density needed for screening the longitudinal electric field. These two parameters 
are connected with total jet power Pjet by equation (|BeskinJ2010J) 


Om : 


JVM 

A I Pa) 


1/2 


( 10 ) 
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Here Pa = m^c 5 /e 2 ~ 10 17 erg/s is the minimum energy loss rate of a central engine 
which can accelerate particles to relativistic energies. 

There are two theoretical models of plasma production in a jet. In the first one, 
pairs are produced by two-photon collisions with photons with sufficient energy pro¬ 
duced by the inner parts of the accretion disk ( Blandford and Znajek|1977f . In this 
case we expect X ~ 10 10 -10 13 , and Michel magnetization parameter Om ~ 10—10 3 . 
The second model of the pair production is a cascade process in non-zero electric 


field in a region with the zero GJ plasma dens ity due to general relativity effects (Be¬ 
skin et al. 1992; Hirotani and Okamoto 1998). This model gives X ~ 10 2 -10 3 , and 
magnetization Om ~ 10 1U -10 13 . So in both scenarios the flow is strongly magnetized 
at its base, but this large difference in the estimates for the magnetization parameter 
Cm leads to two different pictures of the flow structure in jets. Indeed, as has been 
shown by several authors (Beskin and Nokhrina 2006; Tchekhovskoy et al. 2009 


Komissarov et al.|2009} , for well-collimated magnetically dominated MHD jets the 
Lorentz factors of the particle bulk motion follows the relation 


T « r ± /R L , 


( 11 ) 


where r is the distance from the jet axis, and R\ = c/Q is the light cylinder radius. 

The relation GD holds until the flow reaches the equipartition regime — the Poynting 
flux is approximately equal to the particle kinetic energy flux. Further acceleration 
is ineffective. For ordinary jets r±/R^ ~ 10 4 -10 5 . As a result, using the universal 
asymptotic solution © one can find that values Om ~ 100 correspond to the sat¬ 
uration regime when there is approximately equipartition between the Poynting flux 
Pem and the particle kinetic energy flux P part . On the other hand, for Om ~ 10 12 the jet 
remains magnetically dominated (P part P e m)- Thus, the determination of the Michel 
magnetization parameter Om is a key point in the analysis of the internal structure of 
relativistic jets. 

The core-shift measurements described above ( |Gould 1979; Marscher 1983; Sokolovsky 
et al.||20lT Pushkarev et al.||2012 1 might provide an observational way to probe X 
and <7m- The ideal MHD flow of plasma with exactly the drift velocity in crossed 
electric and magnetic fields does not emit synchrotron radiation. However, some in¬ 
ternal dissipative processes, for example internal shocks and magnetic reconnection, 
can produce particles with a non-thermal spectrum that can account for the observed 
emission. By measuring the core-shift discussed in the previous subsection, one can 
obtain the combination of magnetic field magnitude B cs and electron number density 
N cs in the region with optical depth equal to unity. On the other hand, B cs and N cs can 


be related through the flow magnetization <7 = Oyj/T — 1 (Tchekhovskoy et al. 2009) 
assuming that the flow is in the saturation regime, i.e. <7 CS « Om/2 ~ r cs (Beskin and 


Nokhrina 2006, Tchekhovskoy et al.|2009 Komissar ov et al.|200 9). Combining the 
latter with the definitions of X and <7m • one can obtain both X and <7m as a function 
of core-shift and total jet power Pjet- The latter can be estimated by correlating the 
radio jet luminosity to total power needed to form the cavities in surrounding gas 
(Cavagnolo et al. 2010). This analysis leads to estimates of X ~ 10 13 and cr ~ 10 
( Nokhrina et al.||2015 ). The measurements of core-shift by (Pushk arev et al.|[2012 ) 
thus allow the order of Michel magnetization at the jet base to be estimated, and it 
can be concluded that the MHD flow on the observable scales are in the saturation 
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Fig. 6 Magnetic field strengths on different size scales illustrated using a few well-studied objects on pc 
and kpc scales. The size scale plotted here is the transverse radius of the component in which the field 
is measured, not the distance from the nucleus (though in general the two scale). Data are taken from 
OISunivan^and^Gabuzdajj^OO^bJJSavolainen^^ik^Wlta^aandTomikoskijj^OOSUHardcastle^Il^ris, 
WorrallandBirkinshaw (2004 1 , Hard castle andCro ston 1 2010), |Croston an dHardcastle 1 20141, |Hardcastle 
and t 'rosIon|j21) ] I ) and Abdo et al.|j2010} . The solid line shows an arbitrarily normalized line of li - x r 
— we would expect the normalization of this line to depend on the jet power so it is not surprising that all 
points do not lie on the line. 


regime, i.e. there is no effective plasma acceleration by the MHD mechanism down¬ 
stream. Moreover, the order of Om correlates well with the estimated Lorentz factors 
of jet flows discussed above. 


Recently, it has become possible to use the core-shift to measure the magnetic 
flux threading the jets and central black holes in a rather large sample of radio-loud 
active galactic nuclei (Zamaninasa b et al.|2014||Zdziarski et al.|2014) . Dynamically- 
important black hole magnetic fields were inferred, indicating that the central black 
holes of most if not all radio-loud AGN are in a MAD state (Section [4.2| >. MADs 
have also been inferred in tidal disruption events (Tchekhovskoy et al. 2014[ > and 
core-collapse gamma-ray bursts (Tchekhovskoy and Gianni os||2014| i. This suggests 
that MADs are perhaps not rare or unusual as their name might imply, but possibly 
quite the opposite. 
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5.7 Observational summary 

To summarize, we now have relatively good information about the magnetic field 
strengths on kpc scales, and excellent information about the vector properties of the 
field, which, however, are probably not strongly connected to the field properties 
at jet generation. On pc scales and below, magnetic field strengths are not as well 
constrained, but such constraints from, e.g., core-shift observations as we have are 
consistent with the expectations from models, while polarization observations are 
certainly consistent with the presence of a helical field structure on small scales. 

The field strength constraints obtained for some of the objects discussed above 
are summarized in Fig. [6] which shows a broad consistency with an r ~ 1 scaling over 
the observable ~ 7 orders of magnitude in physical scale. Fields of the order 10 4 
G, and thus magnetic fluxes ~ 10 33 G cm 2 would be implied by an extrapolation to 
scales comparable to the gravitational radius of the central black hole. In principle, 
magnetic fluxes at this level could be dynamically important (plasma j3 ~ 1) and so 
could affect the flow of material in the disk as discussed in Section l4~2l 


6 Conclusions 

Accretion disks and the jets they produce are fundamentally magnetic systems. A 
great deal of progress has been made in the theory of these systems, and in obtaining 
their properties from observation. We began by considering how the entire complex 
system is broken into different components for theoretical analysis. Where then do 
we stand with those components? 

First, accretion disks are turbulent, and that turbulence is almost certainly created 
by the instability to the MRI. This MRI-induced turbulence produces internal stress 
that can transport angular momentum at the rates required. 

Although it seems clear that disks are magnetic, jets appear to require a more 
organized magnetic field than is produced by MRI turbulence. The origin of the jet 
magnetic field remains uncertain. A disk dynamo remains a possibility, and dynamo 
action has been seen in some simulations, but the field produced is not the type re¬ 
quired for jet launching. In the case of stellar jets, the central star may also be the 
source of the required field through its own internal dynamo. The possiblity remains 
that the required field could be carried in from large radii; however, MRI turbulence 
does not seem to accomplish that on its own. More direct infall, taking place outside 
the turbulent disk proper, or in a hot, thick accretion flow with relatively large infall 
velocities, could work if the resulting field is able to reconnect as needed to estab¬ 
lish the required topology on the central star or black hole, or through the central 
accretion disk. 

The general model of jets as magnetized, accelerated and collimated disk winds 
is well established today by observations and theoretical investigations. While for 
young stars the disk wind seems to play the major role in jet launching (Blandford- 
Payne), the jets from AGN and other relativistic sources are plausibly dominated 
by processes energetically supported directly by the spinning black hole (Blandford- 
Znajek, |Blandford and Znajek|1977j l, and observations, at least on the smallest scales. 





30 


John F. Hawley et al. 


show good consistency with the expectations of this type of model, as discussed in 
Section [5] There are still many unknowns. For example, the matter content of rel¬ 
ativistic jets is as yet unclear, and in many cases their propagation speed is poorly 
constrained; there is no direct measurement of the magnetic field of jets from stars or 
of the accretion disk. 

Global MHD simulations are able to treat the global jet launching that is the 
accretion-ejection process. Usually an initial large-scale magnetic field has been as¬ 
sumed so far in these models, only a few simulations were able to consider mean-field 
disk dynamo. How such global, axisymmetric mean-field models match the disk mi¬ 
crophysics and the respective local disk simulations, is still an open question. The 
magnetized disks considered in jet launching simulations use a rather simple model 
for the (turbulent) magnetic diffusivity or viscosity and also neglect radiation and ra¬ 
diative transport processes. However, all these processes seem to be essential for the 
disk physics. 

Future progress will clearly come from more sophisticated simulation models, 
considering the physical effects mentioned above. The codes and computational re¬ 
sources may not yet be ready for such an ambitious goal, as the treatment of the 
additional physics will require new numerical methods, and also a higher numerical 
resolution, together with a fully 3D treatment. 

Simulations have clearly demonstrated the viability of the Blandford-Znajek type 
jet powered by a rotating black hole. Those simulations further suggest that a dipole 
field is required, in addition to the rotating hole, and that the field strength is set by 
equipartition with the surrounding pressure (gas, ram, magnetic) associated with the 
accretion disk. The issue of the presence of such jets comes back to the origin of the 
axial black hole field, a problem that remains uncertain, as discussed above. 

Finally, we have seen that observations, particularly of AGN jets, are starting to 
provide direct tests of jet-launching models through constraints on the bulk speed, 
particle content, kinetic power, and magnetic field strength and configuration of jets. 
Much progress on at least some of these issues can be expected from the improvement 
in observational capabilities at arcsecond resolutions (and thus generally kpc scales) 
to be provided by next-generation radio facilities such as the upgraded JVLA and 
LOFAR, and in the future the SKA. However, the most direct constraints remain 
those provided by long-baseline interferometers such as the VLBA and EVN, and 
continued support of these is vital if we are to be able to make direct tests of jet 
models in the future. 
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